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The Isometric Fluctuation Relation (IFR) [P.I. Hurtado et al., PNAS 108 , 7704 (2011)] relates 
the relative probability of current fluctuations of fixed magnitude in different spatial directions. 
We test its validity in an experiment on a tapered rod, rendered motile by vertical vibration and 
immersed in a sea of spherical beads. We analyse the statistics of the velocity vector of the rod and 
show that they depart significantly from the IFR of Hurtado et al. Aided by a Langevin-equation 
model we show that our measurements are largely described by an anisotropic generalization of 
the IFR [R. Villavicencio et al., EPL 105 , 30009 (2014)], with no fitting parameters, but with a 
discrepancy in the prefactor whose origin may lie in the detailed statistics of the microscopic noise. 
The experimentally determined Large-Deviation Function of the velocity vector has a kink on a 
curve in the plane. 

PACS numbers: 45.70. -n, 05.40.-a, 05.70.Ln, 45.70.Vn 


Fluctuation relations (FR) go beyond the second law 
of thermodynamics by quantifying the relative probabil¬ 
ity of short-time, small-scale entropy-consuming events 
|l|. Experimental evidence for their validity is surpris¬ 
ingly widespread even in non-thermal noisy sys¬ 

tems nominally outside their realm of applicability [61-4111. 
A more stringent symmetry property of thermal systems, 
the Isometric Fluctuation Relation (IFR) [l2[, has been 
derived recently for isotropic systems, describing the rel¬ 
ative probabilities of observing currents of equal magni¬ 
tude in different directions, not necessarily diametrically 
opposed as in the standard FR [l|. 

In this Rapid Communication we present the first ex¬ 
perimental observations of behaviour consistent with the 
IFR, on a macroscopic, fore-aft asymmetric rod which 
executes self-propelled 13|-18 motion through a back¬ 
ground of non-motile spheres. Our main results are as 
follows: (i) Our anisotropic experimental system de¬ 
viates substantially from the predictions of [12]. (ii) 
We show that the symmetry properties of the Large- 
Deviation (LD) function correspond to the Anisotropic 
IFR of 0 , with which we are able to show a parameter- 
free agreement through a simplified description based on 
an ansotropic single-particle Langevin equation. There 
remains a 20 % discrepancy in the prefactor whose ori¬ 
gin, we speculate, could lie in non-Gaussianity in the 
microscopic noise, not accounted for in our Langevin 
model, (iii) We find that the measured LD Function 
of the velocity vector exhibits a kink where the velocity 
component along the rod axis vanishes. We speculate 
why these behaviours should arise in an apparently non¬ 
time-reversible system. 

In our experimental setup [lf| a single geometrically 
polar brass rod, 4.5 mm long with diameter 1.1 mm at 
the thick end is placed amidst a monolayer of spherical 


aluminium beads of diameter 0.8 mm (Fig. [H a ))- The 
beads lie on an circular aluminium plate 13 cm in diam¬ 
eter, covered by a glass lid at 1.2 mm above the surface, 
thus forming a confined two-dimensional system. The 
bead area fraction, based on their projected images, is 
0.83. A permanent-magnet shaker (LDS V406-PA 100E) 
drives the plate sinusoidally in the vertical direction with 
amplitude ao and frequency / = 200 Hz, corresponding 
to dimensionless shaking strength T = ao(2Trf) 2 /g = 6.5, 
where g is the acceleration due to gravity. The rod trans¬ 
duces the vibration into predominantly forward motion 
as indicated by the arrow in Fig. UTa) 0Q0. The 
bead medium is both athermal noise source and obstacle 
course for the motion of the polar rod (see Supplemen¬ 
tary Video). A high-speed camera (Redlake MotionPro 
X3) records the dynamics of the particle at a rate of 50 
frames per second. A typical experimental snapshot is 
shown in Fig. [lja) . The images were analysed in Im- 
ageJ 0 to calculate the instantaneous in-plane position 
R(f) and orientation, i.e., the unit vector ri||(f) from the 
thick to the thin end of the rod of the polar rod, and its 
two-dimensional velocity vector v(f) at time t defined as 
the discrete time derivative of R (t) between successive 
frames. Note that v(f) is already coarse-grained in time 
with respect to the true microscopic velocity. The plate 
vibrates at 4 times the frame rate, and the collisions of 
the rod with plate, lid and beads take place at irregular 
instants. 

We resolve v(f) into components (i>|| (f), v±(t)) = (v(f)- 
n||(f),v(t) • nj_(t)) along and transverse to rf||(f) [Fig. 
[He)]. In earlier work [16[ on this system we studied the 
statistics of U|| alone and found a large-deviation func¬ 
tion with a kink at zero as predicted in several theoret¬ 
ical treatments of forced Brownian motion amidst pe¬ 
riodic obstacles [2II-I24I]. The results of 0 encourage 
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FIG. 1: (a) Geometry of the polar particle: the thick 
arrow along the particle’s axis indicates the mean 
direction of “self-propelled” motion and typical 
experimental screen-shot depicting the particle moving 
its way through the crystalline medium consisting of 0.8 
mm aluminium beads in a confined two-dimensional 
cell, (b) Time-series of velocity components U|| & v± of 
the particle. Positive direction is along the green arrow 
in (a). Probability distributions of t>||, (c) and v ±, (d) 
with solid lines showing a Gaussian fit. (e) Resulting 
probability distribution of velocity v = uiiri|| + w_i_rij_ of 
the particle. Circle shows a constant-velocity contour 
and the arrows representing isometric vectors 
subtending angles 9 & O' with the horizontal axis. 


us to further probe its nonequilibrium fluctuations and 
look for a possible IFR. Fig. [IJb) shows the instan¬ 
taneous time-series of U||, whose nonzero mean signals 
systematic “self-propulsion”, and v±, with mean zero. 
The statistical anisotropy of the dynamics is clear from 
the probability distributions of U|| and v± in Figs. [TJ c) 
and (d) respectively, which show a much greater disper¬ 
sion along the rod than transverse to it [25[. Fig. 0e) 
shows the distribution P(v) of the two-dimensional ve¬ 
locity vector, peaked at a non-zero U|| with significant 
weight in all directions including backwards. Our exper¬ 
iments can thus explore the applicability of the IFR in a 
much larger angular range than in the numerical study of 
0 - Before presenting our experimental findings in de¬ 
tail, we build a minimal single-particle Langevin equation 
model for the dynamics of the two-dimensional position 


R(f) = (X(t),Y(t)) of the polar rod as a function of 
time £, ignoring inertia. This will allow an independent 
determination of parameters required later in the paper, 
and help relate our findings to the anisotropic [19| iso¬ 
metric fluctuation relation in a simplifying limit. As we 
work in a frame fixed in the particle, with X and Y de¬ 
fined with respect to the particle orientation, the dynam¬ 
ics of ri|| (f) does not enter our analysis. The governing 
Langevin equations read 

r R = F + N f (t) (1) 


where F, T are the systematic force and damping matrix. 
F gets contributions from the propulsive force driving the 
particle along ri|| (t) as well as from the spatial structure 
of the obstruction to motion posed by the bead medium. 
The random kicks that the rod receives from beads and 
plate are described by a noise f. The physical noise co- 
variance is N-N = 2DTT, which defines the diffusion ten¬ 
sor D. We assume the unit-strength noise f to be Gaus¬ 
sian and white, a plausible assumption whose validity can 
be tested only by comparison to experiment. Our interest 
here, as in [12j . is in large deviations of the macroscopic 
current J J (r, t ). For the present single-particle system 
the current density J(r, t) = d(r — R(f))v(t) at point r 
at time t , where v = R, so that the macroscopic current 
is simply v(f). 

We would like to extract the parameters in 0 from 
the experiment. We define V r (f) = r _1 f* +T v(t')dt', the 
current coarse-grained on a timescale r. For the smallest 
accessible r, of order the inverse frame rate, it is reason¬ 
able to suppose that the structure of the bead medium 
does not affect the dynamics of the rod significantly. We 
can then assume F in 0 depends only on time. Define 
n = N -1 r, S = N^F, and S T (t) = r” 1 / t * +T S(t')dt', 
we can construct, from 0 and the statistics of the noise 
sources, the probability density 


Pr(Vr 


/ T \ d / 2 

V) = detn (-) exp 


--(n- v 

. 2 V 



for V r (£) to take a value V = (Vj|, Vj_), via its moment¬ 
generating function as shown in Appendix. If we ap¬ 
proximate S T (t) by n v 0 , where v 0 is the steady-state 
average velocity, which would be exact for r —>■ oo and is 
reasonable for the present case where the inverse of frame 
rate, r/, is four times the oscillation period, 0 becomes 


P T (V T = V) = (detD)~ 1/2 (-^ d/2 e -i( v ~ v o) TD I ( v - V °) 

(3) 

where the inverse diffusion tensor (2D) 1 = n T II can 
be seen to provide a natural inner product. Below we 
use the form 0 to extract values for uo and the diffusiv- 
ities from our data. Eq. 0, the result of ignoring the 
position dependence of the forcing in 0 , trivially obeys 
an anisotropic IFR 13 because its large-deviation func¬ 
tion, from 0, is quadratic: two coarse-grained currents 
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V and V' satisfying 

V t D' 1 V = V ,t D _ 1 V' (4) 

i.e., which lie on the ellipse V T D _1 V = constant, obey 
1,. ^(V T =V) 


™ T ln P r (V r = V') 


= e .(V-V / ). 


(5) 




with 


e = v 0 r (2D)- 1 . ( 6 ) 

In the event that DT turns out to be proportional to the 
unit tensor, with coefficient T e ff , local detailed balance 
holds. Then e = rv 0 /2T e // is the drag force scaled by 
effective temperature, and (0 becomes a true fluctuation 
relation for the power. 

For comparison with our experiments, let us consider 
currents V, V' with equal magnitude V. The result of 
[Til ] can then be re-expressed, for the case of a diagonal 
D = diag(L>||,£)j_), as 


1 , Pt(V t = V) 

; ln p T (v, = vo = 


2D\ 


■(cos 9 — cosd') 


Vn 


2 D± 


(sin 9 — sin O') 


V 2 (D\l-D±) 

4-D_l-D|| 


(cos 2 9 — cos 2 9') 


(7) 


Although the large-deviation function for our system will 
not have the quadratic form implied by 0 , the foregoing 
calculation gives us a value [Eq. 0] for e in 0 in terms 
of independently measurable quantities. This allows a 
parameter-free comparison of our measurements to the 
Anisotropic Fluctuation Relation (AIFR) 0 in the form 
©or 0 . 

We return now to the experiments. The polar rod 
is propelled only along its nose, i.e., in the x direc¬ 
tion. Thus y-component of F in 0 is zero, so that 
only the component of the mean velocity is nonzero. 
From the measured time series of positions we can ex¬ 
tract the distribution of V T for the shortest time accessi¬ 
ble (77) and infer i>q, D\\ and D± by fitting the mea¬ 
sured distribution to ©. We obtain D\\ and D± as 
((V||(t) --Co) 2 ) t T)/2 and (({V±(t)) 2 ) t r)/2 respectively 
where r = 17 = 0.02 s and vq = 0.34 cm s -1 . We now 
plot constant velocity contours centered at (Vj| = 0, V± 
= 0) and consider sets of V T and V(. at angles 9 and 
9’ with respect to rf|| (Fig. 0e)). We consider overlap¬ 
ping azimuthal bins of V r in order to improve statis¬ 
tics, and obtain the probability density P(V T ). Defining 
J V(D l{ - D ± ) 

0 ,d') = + Jl ±; (cosfl + cosd'), we plot 

2D || 4D_lD|| 

(flr) _ 1 ln[P T (V)/P r (V')] as a function of cosd - cosd', 



FIG. 2 : (a) A typical plot of (Hr )" 1 ln [P T (V)/P T (V')] 
vs. cosd - cos O' over various constant-velocity contours 
for r = 0.26 s showing a linear trend for all V. (b) Data 
scaling of (HrF ) -1 I 11 [P r (V)/P T (V')] vs. cosd - cosd'. 

(c) Scaling of (DrF ) -1 ln [P T (V)/P r (V')] with r 
variation. Here each r line contains all the V values as 
in (b). (d) Plot of (2Drcosd )- 1 ln [P T (V)/P r (V')] vs V 
for various r for the special case when 9 — 9' = 180°. 
Here d varies between —30° to 30° in the step of 10° for 

all r. 


as shown in Fig. 0a) for a typical r = 0.26 s at vari¬ 
ous values of V. The trend is linear for all V in agree¬ 
ment with Eq. 0 (please recall that =0). Further¬ 
more, a plot of (HrI/) _ 1 ln[P r (V)/P r (V')] shows a clear 
scaling with V for r = 0.26 s in Fig. 0b). Similar 
data collapse is seen for all r > 0.10 s. Now we plot 

for various r in Fig. 0c). 
Here each r line contains various V values as in Fig. 
0b). Overlapping lines for all r confirm the r scaling 
of the AIFR. From Eq. 0 the expected slope in Fig. 
0b) or (c) shold be unity. Our experiment, using the 
estimates for vq, D\\ and D±, finds a number close to 
0.8. Perhaps this reflects limitations inherent in using 
estimates of mean speed and diffusivity from short-time 
data. A possible reason for this discrepancy is that the 
noise is not Gaussian [2^]: Fig. 0d) shows significant 
non-quadratic tails in the logarithm of the distribution 
of the small-r averaged velocity, giving some credence to 
this suggestion. 

As a check, we consider the special case of op¬ 
positely directed vectors, i.e. V'=—V, (implying 
9 - 9' = 180°). Fig. 0d) shows the plot of 
(2Drcosd) _ 1 ln[P r (V)/P r (V')] as a function of V for var¬ 
ious r. Here d 6 (—30°, 30°) separated by a step of 10°. 
A clean collapse of data is observed with the same slope 
as in Fig. 0c). This illustrates compliance with the 
standard FR [l|] as a special case of AIFR. 
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FIG. 3: (a) Ignoring the anistropy term in Eq. [71 i.e. 
assuming D\\ = D±, no scaling is observed for (tH) _1 
In [P T (V)/P T (V')] vs. cos 9 - cos 9'. (b) AIFR analysis 
for the particle moving on a bare substrate in the 
absence of bead medium. No collapse of data for 
(VLtV)~ 1 In [P r (V)/P r (V')] vs. cosd - cos 9' suggesting 
the presence of the noisy medium for AIFR to hold. 


For completeness, we analyse the velocity fluctuations 
in the experiment in terms of the IFR for isotropic sys¬ 
tems [12]. It can be obtained by setting D\\=D± in Eq. 


[7] which gives tt = Vo/2D\\ independent of 9,6'. We plot 
(ry)” 1 ln[P r (V)/P r (V')] against cos 9 - cos 9' for r = 
0.26 s. Fig. (3[a) shows that the isotropic IFR is not sat¬ 
isfied. We check that this trend is seen for all r. We 
find a trend of increasing slope with V which cannot 
be explained without the inclusion of the second term 
in fl. This reiterates, experimentally, the importance of 
anisotropy as pointed out by 0 . We check in addition 
the role of the noise and hindrance to motion provided by 
the bead medium, by examining the velocity fluctuations 
of the particle for its motion on a bare substrate in the 
absence of any medium. Now the source of all the noise is 
multiple random collisions of the particle with the base 
and the lid. After an analysis along the same lines as 
above, we plot in Fig. E[b) (GTy)- 1 ln[P r (V)/P r (V')] 
as function of cos 9 - cos 9' for r = 0.12s. Interestingly, 
we find a significant failure of data-collapse in this case. 
Similar behaviour is observed for all r. The presence of 
the bead medium appears to be important for the AIFR 
to hold, a result consistent with our earlier findings [l 6 | . 


It is curious that fluctuation relations, which are 
a derived consequence of microscopic time-reversibility, 
should arise in a system with a unidirectional flow of en¬ 
ergy. It would appear that in the present experiment the 
only significant effect of this energy input is its transduc¬ 
tion into directed motion by the polar rod. The central 
role of the bead medium is also a puzzle. We speculate 
the medium provides closely related processes governing 
the damping and diffusion matrices T and D in (|T|) of the 
polar rod, thus giving rise to an effective local detailed 
balance. In addition, perhaps the multiple collisions of 
the polar particle with the bead medium provides sup¬ 
press correlated movements of the rod due to rolling or 
sliding, which may be present in motion on a bare sub¬ 
strate. 



V„ (cm/s) V x (cm/s) 


FIG. 4: (a) The Large Deviation Function for the 
velocity fluctuations for the case of the polar rod 
moving through the bead medium is non-paraboloid, 
(b) The side-view along the Vj_-axis shows an abrupt 
fall which resembles the kink-like feature observed in 
0. (c) The side-view along the Vj|-axis. 


Lastly, we extract the Large Deviation Function (LDF) 
for the velocity vector of a polar rod moving through the 
bead medium, G(V) defined as P T (\p = A T exp(rG(V)) 
(G(V) < 0) in the limit r —> oo [27]. Here we take 
A r = max[P T (V)] = P T ((V) = vq) which is indepen¬ 
dent of V as required. We plot the asymptotic LDF as 
shown in Fig. |U]a). Clearly the shape of LDF is far from 
paraboloid. There is a sharp drop in G(V) below Vj| = 0 
plane which is clearer when viewed along the Vj_-axis 
(Fig. [U]b)). This is related to the kink observed in ear¬ 
lier work fl 6 |. Fig. me) shows that the LDF is parabolic 
when viewed along the V± axis. 


We point out here that earlier tests of the IFR, e.g., 
energy diffusion on a two-dimensional lattice, a hard¬ 
disk fluid in a temperature gradient [ 12 ] were numeri¬ 
cal. Moreover, these simulations saw no negative events 
and fluctuations were only along the imposed gradient, 
resulting in a test of the IFR in a very limited range of 
angles 9. Our experimental study, with substantial noise, 
allows an exploration over a large 9 range. 


To summarise: we find detailed experimental support 
for the anisotropic variant 0 of the Isometric Fluctu¬ 
ation Relation [12[, for a self-propelled granular particle 
moving through a monolayer of spherical beads. We sug¬ 
gest that the discrepancy between the observed and pre¬ 
dicted prefactor arises from non-Gaussianity in the form 
of heavy tails in the microscopic noise, not accounted for 
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in our Langevin equation model. More work is needed to 
address this. Our measurements are consistent with an 
earlier study of a Fluctuation Relation, and include an 
extension of our earlier observation of a kink in the Large- 
Deviation Function of the velocity vector of the particle. 
Our treatment through an anisotropic Langevin equation 
sheds some light on the circumstances in which such re¬ 
lations should arise, and allows a parameter-free test of 
the theory. 

For support, NK thanks the University Grants Com¬ 
mission, India, HS thanks the CSIR, India, and AKS and 
SR acknowledge a J C Bose Fellowship of the DST, India. 


APPENDIX: Detailed calculation for Anisotropic 
isometric Fluctuation Relation 

The overdamped Langevin equations for a particle sub¬ 
jected to a force F (t) in d dimensions can be written as 

r-R = F + N-f(t). (8) 

Multiplying above Eq. by N 1 gives 


II-R = S + f(i), (9) 

where II = N _1 T, S = N X F. The current for the single 
particle at point r will be 

J(r, t) = S(r — R(t))v(f). (10) 


The macroscopic current averaged over time r then 
reads 


v r (t) = T 1 J t ^J(M) 

= r ” 1 J J S(r — R(t'))v(t')dt' 


= r 


-l 


rt+r 

J v{t’)dt'. 


( 11 ) 


Average of the Eq. © over time duration r gives 


n-v T (t) = s r (i) + J>(i), ( 12 ) 

where S T (t) = r _1 / t * +r S(t’)dt' and -F r (f) = 

r ” 1 / f t+r i(t')dt'. The probability density for F T (t) can 
be written as 


P{Tr(t) = A) = ( S 

i / r (i 

exp —i 




(2tt) 


rt+r 


k • f (t')dt' - k • A 


(27r) fl 


J ^exp —i—Jt k f {t')dt^j expik-A(13) 


But from a well known identity 


exp I / H(t)g(t)dt 


= exp 


H{t) 2 dt 


(14) 

for any arbitrary function H (f) and any white Gaussian 
noise g(t) with zero mean and deviation one, 


P{Tr(t)= A) = 


1 


(27r) fl 

1 

(27r) fl 


exp 

exp l - 


k 2 


k 2 dt + *k • A 


2 r 


Ik • A 


(15) 


RHS of above equation is just Inverse Fourier Transform 
of a Gaussian function exp — fc 2 /2r, therefore the proba¬ 
bility density for !F T {t) becomes 


JWO-A,-(£)*■.,(-If). 


tA 2 


(16) 


From Jacobian transformation the probability density for 
V r (t) is 

'dF T 


P(V T (t)) = det 


dV 7 


P(Tr{t)). 


But from Eq. oa 


dTr 

Wr 


n. 


(17) 


(18) 


hence 


P(V r (f) = V)=detn(^-) 


d/2 


exp 

Since n T II = (2D) 1 , 


(_l (n .V-5 T (t)) 2 ). (19) 


and 


detn = (det(2D~ 1 )) 1 / 2 
= 2 _d / 2 (detD) _1/2 

(II-V-<S r (f)) 2 

= (n v - s T (t)) T (n v - s T (t)) 

= (v - rr 1 • s T (t)) T n T n (v - n _1 5 r (t)) 

= I (v - n- 1 • s T (t)) r d 1 (v - n~ 1 5 r (t))(.20) 


Eq. (1T51) then becomes 

PrCVr = V) = (detD )- 1/2 (^) d/2 X 
x(v-n- 1 .s T (t)) r D- 1 (v^n- 1 -s T (t)) 


( 21 ) 


If the average velocity of the particle in steady state is Vo, 
the S T (t) can be approximated by II-vo, thus probability 
density for V T (f) becomes 

Pt (V r = V) = (detD )- 1/2 (^) d/2 X 
e --(V-v 0 ) T D- 1 (V-v 0 ) 


( 22 ) 
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In 2 dimensions 


In two dimensions, for diagonal D = diag(Z?||, D±), Eq. 
(l22l) becomes 


PriYr 


V) 


exp 


T 

47r \J D\\D± 


T 


(V^o) 2 

4^11 


r(V± -^p-) 2 
4 D ± 


(23) 


Two coarse-grained currents V and V' satisfying the con¬ 
dition 


E± 2 

< 

vf 

D ± 

hi 

D± 


(24) 


i.e., which lie on the ellipse V\\ 2 /D\\ + Vj_ 2 /D± = 
constant, obey 

lim - In = V) = £ . (v _ V '). (25) 

r—>oo T P T (V T = V ) 

For currents V, V' with equal magnitude V the result of 
0 can be re-expressed as 
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lim — In 

T—y OO T 


Pr(y T 

Pr(Vr 


11 

V') 


= V 


2D 


-(cos 9 — cos 9') 


2D ± 


(sin# 


+ 


V 2 (D\\ ~ D±) 

4D±D\\ 


(cos 2 9 — cos 2 9') 



(26) 


The true large-deviation function will presumably not 
have the quadratic form implied by (03l) . One should 
note that this calculation gives us a value 

e = (vl/2D h v^/2D ± ) (27) 


for e in (1251) in terms of independently measurable quan¬ 
tities. 
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